Full counting statistics of renormalized dynamics in open quantum transport system 
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The internal dynamics of a double quantum dot system is renormalized due to coupling respec- 
tively with transport electrodes and a dissipative heat bath. Their essential differences are identified 
unambiguously in the context of full counting statistics. The electrode coupling caused level de- 
tuning renormalization gives rise to a fast-to-slow transport mechanism, which is not resolved at 
all in the average current, but revealed uniquely by pronounced super-Poissonian shot noise and 
skewness. The heat bath coupling introduces an interdot coupling renormalization, which results in 
asymmetric Fano factor and an intriguing change of line shape in the skewness. 
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I. INTRODUCTION 

Electronic dynamics of an open quantum system has 
vital roles to play in the development of nanoelectron- 
ics [IHl]- Tremendous effort has been invested to char- 
acterize and control coherent dynamics in a variety of 
nanoscale transport systems A lot of these stud- 

ies are devoted to the energy relaxation and decoher- 
ence resulted from the coupling with a noisy environment 
10 12]. Yet, another important consequence of system- 
environment coupling, i.e., internal energy renormaliza- 
tion, is attracting rising attention recently, owing to its 
essential influence on the internal dynamics of various 
nanoscale systems, such as spin valves |13r U6[ . quantum 
dot (QD) Aharonov-Bohm interferometers [13, Ha], dou- 
ble quantum dots fl9l |20| , and the quantum state under 
measurement [U |22| . 

Full counting statistics (FCS), which characterizes the 
correlations between charge transport events of all orders 
[23|,[24]], serves as a 'spectroscopic' tool to reveal various 
internal mechanisms. In particular, owing to the develop- 
ment of highly sensitive on-chip detection of individual- 
electron tunneling technique, all statistical cumulants of 
the number of transferred particles can now be extracted 
experimentally. Electron transfer statistics in various 
nanoscale systems, such as QD (25rl29j . gold atoms [30j . 
and molecules jU HH , has been investigated and utilized 
to characterize the electronic dynamics. This offers new 
opportunities to study many-body effects such as the 
Coulomb blockade and Kondo dynamics with essential 
implications to nanoelectronics. 

The main purpose of the present work is to relate the 
higher cumulants to the renormalized dynamical prop- 
erties that cannot be accessed by the average current 
measurements. To have a specific example, we will inves- 
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tigate the energy renormalizations in electron transport 
through a double quantum dot system (33|, where QDl 
is directly tunnel-coupled to the electrodes while QD2 is 
side connected to QDl, as schematically shown in Fig.Q] 
The double dots are further coupled to an inevitable dis- 
sipative phonon environment (not shown explicitly) . The 
system is of particular interest, as it is in such a configura- 
tion that maximizes locality versus nonlocality contrast, 
and can be mapped onto nanostrucures in experiments 
011. 

In context of FCS, we are capable of identifying the 
essential difference between the internal energy renor- 
malizations induced by the coupling to the electrodes 
and that due to external phonon bath. The former 
one gives rise to a renormalized level detuning, which 
enhances both Fano factor and skewness uniquely to a 
super-Poissonian value, whereas does not affect the av- 
erage current at all. For the latter, i.e., the phonon bath 
coupling, it introduces an intriguing interdot coupling 
renormalization, causing strong electron localization in 
QD2. Furthermore, the asymmetry in phonon absorption 
and emission leads to an asymmetric shot noise spectrum 
as well as a change of line shape in the skewness. 

The remainder of this Letter is organized as follows. 
Based on the number-resolved reduced density matrix, 
the theory of FCS is outlined in SectionlLTl We then 
introduce in Section llLTl the model Hamiltonian for the 
transport double dots coupled to a heat bath. The phys- 
ical origin of the internal energy renormalizations is dis- 
cussed in Section|lVl along with its influence on the FCS 
analyzed. Finally, all the results are summarized in Sec- 
tionEl 



II. FCS FORMALISM 

Full information about transport properties of a 
given system is contained in the probability distribution 
P(N,t c ) that N electrons have transferred through the 
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system during the counting time t c . The FCS yields all 
the cumulants (/ fc ) of the probability distribution, which 
are directly related to important properties of the junc- 
tion. For instance, the first cumulant (I) gives the av- 
erage current. The shot noise is related to the second 
cumulant (I 2 ) and is commonly represented by the Fano 
factor F = (J 2 ) / (I) , with F < 1 indicating a sub-Poisson 
noise, F = 1 a Poisson fluctuation, and F > 1 a super- 
Poisson process. The third cumulant (I 3 ) measures the 
skewness of the distribution. 

The cumulants can be found by performing derivatives 
of the cumulant generating function (CGF) T(x) with 
respect to the counting field \i 



(I k ) = --j-(~id x ) k F( x ,t c )\ x=0l 



where the CGF is defined as 



^P(N,Qe 



(1) 



(2) 
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The involving probability distribution is closely related 
to the dynamics of the reduced system p(t), which is 
obtained from the density matrix of the entire system 
by integrating out the reservoir degrees of freedom. To 
characterize the FCS, we further unravel p(t) into compo- 
nents p( N \t), in which W denotes the number of elec- 
trons transferred through the system. The probability 
distribution is obtained via P(N,t) = TrpW(i), where 
Tr(- • • ) denotes the trace over the reduced system degrees 
of freedom. 

To proceed, instead of using p( N ^(t) directly, we 
shall employ its x _s P ace counterpart, i.e., g(x,t) = 
^2 N p( N '(t)e lNx . The corresponding quantum master 
equation reads formally |36l439j 



d_ 

di 



g( X , t) =-{iC + K x )q( X , i) = C x g( X , t), (3) 



where £(•••) = [H sys , (•••)] is the reduced system Liou- 
villian, 1Z X is the dissipation superoperator, and L x = 
— (i£ + lZ x ), which will be specified later in the localized 
state representation [see Eq. ©]. 

With the knowledge of g(\, t), the CGF is determined 
straightforwardly as F{x,tc) = — ln{Trg(x, t c )}. In the 
zero-frequency limit, i.e. the counting time t c is much 
longer than the time of tunneling through the system, 
the CGF is simplified to (13-S3 



(4) 



F(X,Q = -A m i n (x)tc, 



where A m ; n (x) is the minimal eigenvalue of C x that sat- 
isfies A min (x -> 0) -> 0. 



III. MODEL DESCRIPTION 

Specifically, let us consider a double quantum dot sys- 
tem, as schematically shown in Fig.[TJ where the QD1 




FIG. 1: Schematic setup for transport through a double quan- 
tum dot system, where QD1 is tunnel-coupled to the left and 
right electrodes, while QD2 is side-connected to the QD1. 
The system is further coupled to an inevitable dissipative 
phonon environment (not shown explicitly). 



is connected to left and right electrodes, while the QD2 
is side-connected to the QD1. The double dots are fur- 
thermore embedded in a dissipative phonon bath. The 
total Hamiltonian consists of the double dot system, the 
environment, and the coupling between them, 



H — H sys + H env 



Hsys— cnv • 



Here the first part models the coupled dots 



(5) 



+ J2e=i,2 Uonefriei + U'hih 2 , (6) 



where Q za = d\ a d Xa - dl a d 2rT , Q xa = d\ a di a + d\ a d\ a , 
™t = J2a ™ttr, and h ta = d\ a d iai with d ia (d\ a ) the elec- 
tron annihilation (creation) operator in the QD1 (1=1) 
or QD2 (£=2) and spin a =f or \.. Each dot consists 
of a single spin-degenerate energy level E 1 / 2 — E ± |e, 
measured relative to the equilibrium chemical potential 
of the electrodes. Electron tunneling between the two 
dots are characterized by the interdot tunneling strength 
O. Double occupation on the same dot costs the intradot 
charging energy U$. Simultaneous occupation of one elec- 
tron in each dot is associated with the interdot charge 
energy U' . Hereafter, the intradot charging energy Uq is 
assumed to be much larger than the bias voltage, such 
that charge states with three or more electrons in the 
double dot are prohibited. 

The environment is composed of the phonon bath, the 
left and right electrodes. The corresponding Hamilto- 
nian reads H env = h p h + h e \. Each of them is modeled 
as a collection of noninteracting particles. The phonon 
bath are modeled with /i p h = (1/2) Y2j ^jiPj ~^~ X J)- The 
electrodes assumes h cl = ^ q =l,r J2ka ^ak^ aka c a k a , ex- 
pressed in terms of electron creation (c^ aka ) and annihila- 
tion (cakcr) operators in the left (a =L) or right (a =R) 
electrode. The electrodes are assumed to be in equilib- 
rium, so that they are characterized by the Fermi distri- 
bution /l/r(w)- An applied bias voltage V is modeled 
by different chemical potentials in the left and right elec- 
trodes /iL/R = ±eV/2. 
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The system-environment coupling can be written as 

sys-cl + -Hsys-phj with 



Hsys-el = X^^ QfcC «fccr rfl<:r + H - C 0i ( 7a ) 
-Hsys-ph = ^ QzaFph- (7b) 

(T 

Here H sys ^ c \ describes electron transfer between QD1 
and the electrodes. iJ sys ^ph models the coupling with 
phonon bath, in which F p h = Y) . XjXj. This term is re- 
sponsible for the dot level energy fluctuations. The tun- 
nel coupling between electrodes and QDf is characterized 
by the coupling strength T a — 2irJ2k \tak\ 2 8(ui — £ a k)- 
In what follows, we consider only spin conserving tun- 
neling processes, and assume wide bands in the reser- 
voirs, which yields energy independent couplings F a . 
The effect of phonon bath on the double dot system is 
characterized by the phonon interaction spectral density, 
J p h(o;) = J^j I -^j 1 2< ^( aj — Hereafter, it is assumed to 



be Ohmic, i.e., J p h(w) 



rjuie 



[J ' 0Jc , where the dimen- 



sionless parameter 77 reflects the strength of dissipation 
and lo c is the Ohmic high energy cutoff. Throughout this 
work, we set unit of % = e = 1 for the Planck constant 
and electron charge, unless stated otherwise. 



IV. LEVEL RENORMALIZATION AND FCS 
ANALYSIS 

In the strong intradot Coulomb blockade regime, dou- 
ble occupation on the same dot is prohibited. The in- 
volving states are restricted to: |0)-both dots empty, 
|ler)-one electron in QD1, |2cr)-one electron in QD2, and 
|lcr2(T')-one electron in each dot, respectively. The quan- 
tum master equation (|3]) in this localized state represen- 
tation reads 

go = -2(r+ + T+)g Q + (T£ + T R e ix )(g n + g u ), (8a) 
g la =in{g\l-gil)-{Tl+T K )g la 

+ (r+ + r+e- i *)ft )) (8b) 

g 2a = ifi(g£ - g&) - 2(f + + t+)g 2a 

+ (rL + r R e [X )(gia2a+gi*2a), (8c) 

gi*2«> = -(f l + f r >i ct2ct < + (f + + f+e- ix )g 2 „', (8d) 
q\1 = i(e + e)g\l + i(fi + tt){g la - g 2a ) 

+ r + gi CT - r_g 2o - - r d £>2^, (8e) 

with spin a e {t:l} an( l o = —o- Here g s — (s\g\s) 
represents the diagonal element of the reduced density 
matrix. The off-diagonal elements g s s , — (s\g\s') de- 
scribes the so-called quantum "coherencies" . The involv- 
ing temperature-dependent tunneling rates are defined 
as T± = T a fi ±} (E ) and f± = T a fL ±] (E + U'), where 
f^\uj) = {l+e^^-i 1 ^}- 1 is related to the Fermi func- 
tion of the electrode a= L or R, with (3 = (fceT) -1 . Here 



we are interested in the regime A < k^T (A = \J e 2 + 4il 2 
being the eigenenergy separation), where the external 
coupling strongly modifies the internal dynamics, and the 
off-diagonal elements of the reduced density matrix have 
vital roles to play [43|, [44| . The level separation is thus 
smeared by the temperature and only the excitation en- 
ergy levels Eq and Eq + U' enter the Fermi functions. 

The coupling to the phonon bath results in relaxation 
and dephasing in the system, with corresponding rates 
given by |45l448j 



r± 

and 



(3A 2 



2 A 



^ph(A) 



A 



• coth 



n 



7 P h = 2tt— Jph(A) coth — 



/3A 



■que 



pA 2 



The total dephasing rate eventually reads 



7 P h- 



(9a) 



(9b) 



(10) 



By observing the equation of motion of the off- 
diagonal matrix element [Eq. ([5e])] . it is found that the 
energy scales of the double dots are renormalized, i.e., 
e — > e + e and f2 — > O + Q. The energy renormalizations 
e and Cl are entailed respectively by coupling with the 
electrodes and phonon bath, 

e=<t>(E ) - 2 ^Eo + V) + HEo + U ), (Ha) 



with 



1 2 H 2ir 



a=L,R 

D{u) = e- i ^B x {y, -1) - e l %B x (y*, -1 



(lib) 

(12a) 
(12b) 



Here, "F and B x (y,z) are respectively digamma function 
and incomplete beta function [4j|, with x = and 

The detuning renormalization e has been neglected in 
previous works [43|, 0, l5(il - [52T ] , where the Fermi energies 
of the electrodes are assumed to be far away from the 
electronic states of the dots. The effect of interdot cou- 
pling renormalization, to our knowledge, has not yet been 
addressed in the context of electron quantum transport. 

To clearly address the interdot coupling renormaliza- 
tion, we shall employ the quantum master equation ([3]), 
in which the dissipative term contains the contribution 
from the phonon bath 



K phg (x,t) = [Q z ,Q z Q{x,t)]+H.c. 



with Q z = J2 a Qz°- and 



Q z = [C(-£)+LD(-£)]Q 2 



(13) 



(14) 
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FIG. 2: (a) level detuning renormalization, (b) average current (/}, (c) Fano factor F = {I 2 ) /(I), and (d) normalized skewness 
S — (I 3 ) /(I) versus the bias voltage for different Fl/Tr ratios. Each time when the chemical potential of the electrode aligns 
with the energy needed for either single (Eo) or double occupation (Eo + U'), e reaches its local extremum. Plot parameters 
are: e=0, r = T L + Tr = 2Q,, E = 10fc B T, U' = 25fc B T, and U = 60fc B T. 



Here, C(— £) = dtC(t)e lCt is spectral function. 
The involving phonon bath correlation function is C (t) = 
(F ph (t)F ph (0)) B , with (•••)= Tr B [(---)M and p B the 
thermal equilibrium state of the phonon bath. The dis- 
persion function D(—C) involved can be evaluated via the 
Kramers-Kronig relation D(—C) = — -V f°° du c r ^ , 
where V denotes the principal value. This term is respon- 
sible for the energy renormalization [l2[ . Let us focus on 
the last term in Eq. (IT4"]) , which can be further expressed 
as 



LD(-£)& 



A7ri^[D(A) + D(-A)]Q 



(3A 2 



(15) 

with Q + = |1)(2| and Q_ = |2)(1|. Apparently, the first 
term accounts for the renormalization of level detuning. 
However, due to the odd function D(uj), this term even- 
tually vanishes. The last two terms represent phonon 
absorption and emission. It is this phonon-induced effec- 
tive coupling that gives rise to a unique interdot coupling 
modification which takes place only in Eq. (|8e| . 

Hereafter, we will show that the both e and have 
important roles to play in the reduced dynamics, and 
consequently give rise to intriguing features in the high 
order cumulants of the current distribution. Consider 
first the situation without electron-phonon interaction 
(rj = 0). The renormalization of interdot coupling fi is 
thus zero. The numerical result for e vs bias voltage is 
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FIG. 3: 3D plot of the interdot coupling renormalization Q 
vs bare level detuning e and temperature fc B T. The plotting 
parameters are Tl = Tr = T/2, and rj = 0.04, u) c = 15T for 
the Ohmic phonon bath spectral density. 



shown in Fig.fJJa) for different Fl/Tr ratios. Once the 
Fermi energy of the electrode is resonant with the energy 
needed for single (E ) or double occupation (E + U'), I 
reaches its local extremum. For a fixed total tunneling 
width r = Tl + Tr, e increases whenever the tunnel- 
coupling asymmetry Fl/Tr increases. 
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Let us now investigate the influence of level detuning 
renormalization on FCS. The first three cumulants of the 
transport current are displayed in Fig.JSfb), (c), and (d), 
respectively, as functions of the bias voltage. Here, the 
results with bare level detuning (e = 0) are also plotted 
for comparison. 

The first cumulant, i.e., average current shows a typical 
step-like structure, as displayed in Fig.[^b). In the low 
bias voltage regime (/xl, A*R < Eq), transport is exponen- 
tially suppressed, as the Fermi energies of both electrodes 
are well below the single level E Q . There are only very 
few thermally activated tunneling events, and the current 
Iq s» 0. Transport through the system becomes energet- 
ically allowed when the Fermi level crosses the discrete 
level, i.e., /xl > Eq > /x R . The current rises to the sec- 
ond plateau. The bias voltage, however, is not sufficient 
to overcome the intradot and interdot charging energies, 
and the coupled dots can accommodate at most one elec- 
tron. This is the so-called double-dot Coulomb blockade 
(DDCB) regime [H,|4^]. The stationary current is ob- 
tained by utilizing Eq. ([1} 



h = 



2r L r R 
4r L + r R ' 



(16) 



where the involving Fermi functions are approximated by 
either one or zero. The effective quadrupling of Tl is a 
combined effect of spin degeneracy and DDCB [37j . 

As the bias further increase and cross interdot charg- 
ing energy, i.e., /xl > Eq + U',Eq > /xr, double occu- 
pancy (with one electron in each dot) is possible. It is 
also referred to as single-dot Coulomb blockade (SDCB) 
regime. The current eventually rises to the third plateau 
as shown in Fig.^b), with its value given by 



h 



2r L r R 
2r L + r R ' 



(17) 



where the effective doubling of Tl arises from the spin 
degeneracy and SDCB [371 ] . 

Noticeably, the current is not sensitive at all to the 
level detuning renormalization e. The reason will be pro- 
vided later. In order to further exploit the effect of the 
energy renormalization, we resort to higher order cumu- 
lants. The numerical result of the Fano factor against 
the bias voltage is displayed in Fig.[2fc). At very small 
bias V <C k^T, it is dominated by thermal noise, which 
is described by the well-known hyperbolic cotangent be- 
havior, and eventually leads to a divergence of the Fano 
factor at V — [H, [54[ ■ As the bias increases but still 
below Eq, there are only thermally activated uncorre- 
cted tunneling events, and Fano factor is Poissonian, 
i.e., Fq = 1. In the DDCB regime, electron transport 
through the bare level Eq is energetically allowed. The 
Fano factor is given by 



F x = l- 



8r L r R 2r^(r R + 4g 2 ) 
(4r L + r R )2 + (4r L + t r )*w ' 



(18) 



where the second term describes the suppression of the 
Fano factor below unity, whereas the third term gives a 



positive contribution. Unambiguously, the Fano factor 
is modulated by the detuning renormalization e. This 
shows the shot noise as a much more sensitive tool to the 
system parameters than the average current. 

The level detuning renormalization gives rise to in- 
triguing features in the noise spectrum. By neglecting 
e, the noise in the DDCB regime is constant against the 
bias. It can reach maximally the Poissonian value (for 
the chosen parameters here Tl + T R = 2f2) in the limit 
of large tunnel-coupling asymmetry, i.e., Tl/Fr ^> 1. 
The detuning renormalization enhance the Fano factor 
in a unique way, i.e., the noise are increased at bias close 
to excitation energies Eq and Eq + U'. Remarkably, a 
pronounced super-Poissonian noise is observed provided 
tunnel couplings are sufficient asymmetric, as shown by 
the dotted curve in Fig.^c). 

To further understand the current and noise features, 
we perform a unitary transformation of the entire Hamil- 
tonian, such that the system Hamiltonian is diagonalized 
to (here spin indices are suppressed for simplicity) 



Ia { \+)(- 



-><-]}, 



(19) 



where A = Ve 2 + 4fi 2 , with the detuning renormaliza- 
tion e being appropriately taken into account. The cou- 
pled dot system can now be mapped onto a parallel two- 
level system. Here, the bonding state |+) and the anti- 
bonding state |— ) are defined respectively as 



= sin-|2)+cos-|l), 

9 9 

• = cos 2 I 2 ) _ sin £ I 1 )' 



(20a) 
(20b) 



where 8 is introduced via sm9 = 20/ A and cos 9 = e/A. 
As a result, the electron tunnel-coupling Hamiltonian 
(l7al) is recast to 



ff sys _ ol = ^{t+ fe 4|0><+|+^4|0)<-|}+H.c, (21) 

ah 

with t~^ k — t a k cos | and t~ k — t a k sin |. Apparently, the 
amplitudes for tunneling through the eigenstates | +) and 
|— ) are effectively modulated by the detuning renormal- 
ization e. 

The stationary current in general can be written as 

1 = J ^TM{/ L M-/ R (u,)}, (22) 

where T(uS), the transmission probability, defines the 
probability for an incoming electron to be transmitted 
through the device. For a parallel two-level system in the 
DDCB regime, the transmission probability is related to 



T 



(I^ 2 + M 2 )I4AI 2 

I^Rfe^Rfc I 2 + 2 I^R/c*Lfc| 2 + 2 |^Rfe*Lfc| 2 



4M 2 + |*Rfe| 5 



(23) 
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where the ^-dependence of the tunneling amplitudes can- 
cel out. Eventually, it leads to the current expression of 
Eq. (|16[) . Physically, suppression of electron tunneling 
through the anti-bonding state due to level renormaliza- 
tion is exactly compensated by that through the bonding 
state. The total transmission probability thus remains 
unchanged, and finally the current turns out to be de- 
tuning renormalization independent. 

The Fano factor, however, is more sensitive to the in- 
ternal dynamics, and reveals much detailed information 
about transport. Particularly, in the presence of large e, 
electron tunneling through the anti-bonding state is sup- 
pressed, while that through bonding state is enhanced. In 
the limit where the Coulomb interactions prevent a dou- 
ble occupancy of system, there is competition between 
the two transport channels. Consequently, the slow flow- 
ing of electrons through the bonding state modulates that 
through the fast state, giving rise to a bunching of tunnel- 
ing events and eventually resulting in super-Poissonian 
noise, as displayed by the dotted curve in Fig.^c). 

Different from the Fano factor, the skewness can be 
either increased or reduced by e. In the case of sym- 
metric tunnel-coupling (Tl = Tr), the level renormal- 
ization leads to a suppression of skewness, see the solid 
curve in Fig.^d). As the coupling asymmetry grows, e 
can strongly enhance the skewness to a super-Poissonian 
value, and pronounced noise peaks are observed at bias 
close to Eq and Eq + U' , as shown by the dotted curve 
in Fig.[^d). The peak width of skewness, however, is 
reduced remarkably, compared with that of Fano factor. 

We now investigate the effect of energy renormalization 
induced by external phonon bath. Let us ficus on the 
DDCB regime, i.e., at most one electron can reside on 
the double dot system, and investigate the effect of Cl 
based on FCS. One can always choose an appropriate 
bias voltage, at which the level detuning renormalization 
is zero, i.e., e = 0, such that its influence is eliminated. 
In Fig. [31 Q is plotted as a function of temperature k^T 
and bare level detuning e. It is observed that the interdot 
coupling renormalization is fully negative; it thus reduces 
the total interdot coupling. It reaches a local minimum 
at e = 0, and grows with increasing e. Q is particularly 
notable at low temperature, but falls off rapidly with 
rising temperature. 

The calculated Fano factor and normalized skewness 
are shown in Fig.UJa) and (b), respectively, as a func- 
tion of the bare level detuning at different temperatures. 
For comparison, the results without fl are also plotted. 
In the absence of the phonon bath (rj = 0), both Fano 
factor and skewness are symmetric around e = (see 
the solid curves in Fig.[4|. With nonzero coupling to the 
heat bath, phonon absorption and emission take place, 
which dominate for e < and e > 0, respectively. At low 
temperatures, the asymmetry in emission and absorp- 
tion leads to the observed asymmetric spectra in Fig. [4] 
The coupling to the phonon bath furthermore causes de- 
phasing between the two dot-levels, which generally sup- 
presses the noise [HI, HH . The reduction of Fano factor 
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FIG. 4: (a) Fano factor and (b) normalized skewness in the 
DDCB regime vs bare level detuning for different tempera- 
tures. The other plotting parameters are the same as those 
in Fig.O 



is particularly notable for e < 0, as shown in Fig.QJa). 
On the other hand, electron-phonon scattering results 
in electron localization in QD2 [48]. In particular, for 
e > and at low temperature (e.g. k^T = T) this mech- 
anism dominates. Electrons thus tend to be transferred 
in bunches, which enhances the noise in comparison with 
the case of r\ = 0, as displayed by the dotted curves in 
Fig.Ha).^ 

Finite £1 (due to its negative contribution) suppresses 
electron tunneling between the two dots, and effectively 
increases electron localization in QD2. Thereby the noise 
is enhanced, particularly at low temperature, where O 
is most prominent [see the dashed curves in Fig.[4ta)]. 
Noise enhancement due to weak interdot coupling were 
also reported in Ref. [H}, where counting statistics 
for electron transport through a double-dot Aharonov- 
Bohm interferometer was investigated. However, there 
the divergence is closely related to a separation of the 
Hilbert space of the double dots into disconnected sub- 
spaces that contain the spin singlet and triplet states for 
double occupancy. 

The effect of phonon bath on the third order cumulant 
is displayed in Fig.[4jb) , where the normalized skewness 
is plotted against the bare detuning e. Whenever the 
double dots are coupled to the phonon bath, the skewness 
is enhanced and becomes asymmetric. Noticeably, above 
the resonance (e > 0) a change of line shape is observed, 
i.e., the skewness increases with growing e, rather than 
decreases as in the case of 77 = 0. This intriguing feature 
shows unambiguously the third order cumulant as a much 
sensitive tool to the phonon bath coupling. The presence 
of finite Q furthermore increases the skewness, especially 
at low temperature [see the dashed curve in Fig.SJb)]. 

With the knowledge of energy renormalizations e and 
f2, we are now in a position to consider their combined 
effect. The Fano factor is plotted as a function of bias 
voltage and temperature in Fig.[S] In the lower half of 
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FIG. 5: Contour plot of the Fano factor versus the bias and 
bare level detuning in the DDCB regime. The plotting pa- 
rameters are, Fl = 2Fr and e — 0. The parameters for the 
phonon bath are the same as those in Fig. [3] 



the bias regime, say V/2 < 20r, the level detuning renor- 
malization e is not strong [see Fig.^a)]. The Fano factor 
thus exhibits sub-Poissonian statistics. Increase of e due 
to rising bias voltage results in noise enhancement, and 
consequently super-Poissonian noises are observed in the 
upper half bias regime. Noise is further increased at low 
temperatures, where the interdot coupling renormaliza- 
tion O effectively suppresses tunneling between the two 
dots and electrons tend to be transferred in bunches more 
readily. 



CONCLUSION 



In summary, we have investigated the full counting 
statistics of electron transport through a side-connected 
double quantum dot system, and paid particular atten- 
tion to the intriguing features arising from the renor- 
malized internal dynamics. It is found that, in com- 
parison with the average current which is not sensitive 
at all to the detuning renormalization, the Fano fac- 
tor and normalized skewness are both enhanced to a 
super-Poissonian value in a unique way. Phonon emis- 
sion and absorption take place when the double dots are 
further coupled to an external heat bath, and give rise 
to an intriguing interdot coupling renormalization. Due 
to negative contribution, it effectively localizes electron 
in QD2, which eventually leads to strong enhancement 
of the noise spectra. The asymmetry of phonon emis- 
sion and absorption processes result in an asymmetric 
Fano factor, and a change of line shape in the skewness. 
Our investigations demonstrate unambiguously the im- 
portance of the internal energy renormalizations, which 
thereby should be properly accounted for in transport 
through various nanoscale devices. 
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